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Abstract 

Stark transitions and induced annihilation of antiproton in the collisions of an- 
tiprotonic helium ions (pHe +2 )„; with He atoms at very low energy (~ 10 K) are 
considered in the framework of quantum coupled-channels method taking into ac- 
count all the states with different I at given n ~ 30, including the annihilating ns 
and np-states. Elastic scattering, Stark transitions and induced annihilation during 
collisions are produced by scalar and dipole terms in the interaction. It is shown that 
the most important contribution to the processes comes from the long-range polar- 
ization interaction. Admixtures of the ns and np-states to the states with higher I 
during collisions induce the effective annihilation cross sections for the initial I up to 
15, but don't affect the Stark cross sections for the initial states nearly to circular 
orbits. Total rates of the Stark transitions from the circular orbits with n = 28 — 32, 
averaged over the thermal motion, are compatible with the recent ASACUSA data. 
Isotope effect as well as the dependence on n are also qualitatively agree with the 
experiment. 



1 Introduction 

The experimental discovery and investigation of the metastable antiprotonic states in he- 
lium have opened a new chapter in the study of antiprotonic atoms (see a comprehensive 
review [1] of the results obtained up until 2001). In addition to the fundamental properties of 
the antiproton, these investigations give interesting insights into the interaction of antipro- 
tonic helium with ordinary atoms and molecules. These data stimulated many theoretical 
papers on the mechanisms of antiprotonic helium formation, on collisional quenching of 
the metastable states, on effects of collisional (density) shift and broadening of El-spectral 
lines, etc. New stage of the study of antiprotonic helium by the ASACUSA collaboration 
at the AD beam (CERN) brings further new data that pose specific problems in the theory 
of collisional effects at low temperature (~ 10 K). Thus for the measurements of hyperfine 
splitting of antiprotonic helium levels by a triple laser-microwave-laser resonance method 
[2] , it is important to estimate the rates of collisional transitions between the HFS sublevels, 
and the density shifts and broadenings of the microwave Ml spectral lines [3j. The first 
observation of cold long-lived antiprotonic helium ions (pRe 2+ ) with lifetimes ~ 100 ns 
against annihilation [4] revealed an unexpected n-dependence and isotope effect for effective 
annihilation rates in circular states (I = n — 1) at n = 28 — 32. Earlier considerations of 
collisional processes of (pHe 2+ ) at higher energy could not explain these data. 
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In this paper we consider Stark transitions and induced annihilation in the collisions of 
antiprotonic helium ion with He atoms at very low energy (~ 10 K), 



( P Re +2 U + He - { + He (St&rk) ' (1) 

I (7r . . . tiX) + He (annihilation). 

We use the quantum close-coupling method with the appropriate interatomic interaction. 
The boundary conditions at R — > oo in the channels corresponding to ns and np-states of 
the antiprotonic ion are modified with regard to the annihilation widths of the states. 



2 Interaction between antiprotonic helium ion and He 
atom 

Let us discuss the interaction between a (pHe +2 ) n i ion and an ordinary He atom at low 
energy (E < E a = 27.2 eV). The relative velocity of the collision is small as compared with 
the electron velocity in the atom, d <C D e , where v e ~ Z e jfV a is a characteristic electron 
velocity and v a = e 2 /h is the atomic unit of velocity. (Hereafter we use the atomic system 
of units h = e = m e = 1, except as otherwise indicated explicitly.) The orbital velocity 
of the antiproton in the states with the principal quantum number n ^> 1 is also small, 
v n — Z/n -C 1. At these conditions the heavy particles (p and nuclei) in the system 
are slow as compared with electrons, therefore the electronic variables can be separated 
out within adiabatic approximation. Then the sum of Coulomb interactions between the 
antiprotonic ion (pHe 2+ ) and He atom is replaced by the effective three-body interaction 
V(R, r) that depends on the relative coordinates R of two colliding subsystems and on the 
inner coordinates r of the antiprotonic ion. 

At a large distance between two subsystems the interaction is reduced to the long-range 
polarization term, V(R, r) — > —a/2R A , where a = 1.383 a.u. is a polarizability of the He 
atom. It is well known that the classical particle in a such potential can fall to the center, 
that is treated as a 'polarization capture' with the effective cross section 

a p {E) = n^a/E. (2) 

Then a mean radius of the effective interaction can be estimated as R ~ (a/E) 1 ^ 4 , i.e., at 
least several atomic units at E < 1 eV. On the other hand, a mean size of the antiprotonic 
helium ion can be estimated as f = n 2 / (fiZ), where [i is the reduced mass of the p — A He + 
system, fi ~ 1467m e for A=4 and 1376m e for A=3. At n = 30 we have f ~ 0.3ao R. 
Therefore the three-body potential V(R, r) can be developed in the multipolar series 

V(R, r) ~ V (R) + (d ■ V R )V (R) + (Q 2 (r) ■ C 2 (Q R )) Vf(R), (3) 

where d = — £ir and Q 2u = — £2?" 2 CW(^r) are dipole and quadruple operators of the ion, 
CW(^) = l22,(f2)^/47r/5. The factors £i = (m^ + 2m p )/(m^ + m p ) and £ 2 = {jn\ ~ 
2m 2 ) I (mA + Tn p ) 2 arise owing to a shift of the center of charge from the center of mass. The 
scalar term Vo(R) in ([3]) can produce only the elastic scattering, whereas the dipole and 
quadruple terms can mix the nl states of the ion and cause the Stark transitions (nl nl') 
of the antiprotonic helium ion. 
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The potential Vq(R) is an adiabatic potential of the interaction between the ionic charge 
(+1) and He atom, that is, in fact, the same as the adiabatic potential of the interaction 
between proton and He atom. The latter, according to the calculations [5], can be approxi- 
mated by a sum of the Morse potential Vm{R) and of the polarization long-range interaction 



V P (R) cut off at an intermediate distance, 

V (R) = V M (R) + V P (R) } (4) 

V M (R) = D {exp[-2/3(i? - R e )} - 2exp[-/3(i2 - R e )]} , (5) 

V P (R) = {1 - exp [- 7 (J2 - R e f] } 6{R - R e ), (6) 

For the calculations we take the following values of the parameters [51 E] Dq = 0.075, 



R e = 1.46, (3 = 1.65, 7 = 0.005 a.u. It should be noted that the Morse potential (jSJ) 
contains a strong repulsion at R < R w = R e — /3 _1 ln2. In order to avoid the loss of 
accuracy at R < R^ in the solving of the coupled-channels equations at very low energy, 
we introduce the infinite wall at R w = 1.0399 a.u. 

3 Coupled-channels equations for the collisions (pHe 2+ ) 
ions with He atoms 

In the previous section we noted that electronic variables in the system (pHe +2 ) n ; + He can 
separated out within the adiabatic approximation. Then total effective hamiltonian of the 
three heavy particles (p — Z — He) is 

H = T{R) + h{r) + V{R,T), (7) 

where T(R) = (— l/2m)V^ is the kinetic energy operator, m is the reduced mass of the 
colliding subsystems, and h(r) is the inner hamiltonian of the antiprotonic ion. For our 
aims it can be presented as 

h = h + U opb (r), (8) 

where ho is the hamiltonian of hydrogen-like atom with the reduced mass \x and with the 
nuclear charge Z — 2, and U opt (r) is a short-range complex optical potential with the 
absorptive imaginary part ImC/ opt (r) < 0. Let E ni and (f> n im(r) = R n i(r)Y [m (Q r ) be the 
eigenvalues and the eigenfunctions of the hamiltonian (jSJ), corresponding to the discrete 
levels of the (pHe 2+ ) system, 

h(f> n i m {r) = E nl (j) nlm (r). (9) 

In addition, for the conjugated hamiltonian h + we introduce eigenvalues E nl = E* nl and 
eigenfunctions n / m (r) = R n i(r)Yi m (Q r ) with R n i(r) = i?* z (r). Two sets of the eigen- 
functions, n «m(r) and 4>nim(r), form jointly a bi-orthogonal system. The eigenvalue E n i 
in general case contains a negative imaginary part, E n i = E^ — ir n ;/2, where r n ; is the 
annihilation width of the level. The real part of the eigenvalue is usually expressed as 
Eni = e n — £nh where e n = —fj,Z 2 /2n 2 is the energy of the hydrogen-like system with the 
point nuclear charge, and (— e n i) is the energy shift of the level due to the nuclear interac- 
tion and the finite distribution of the nuclear charge. Total difference (E n i — e n ) is usually 
referred to as the complex energy shift AE n i = —e n i — iT n i/2. The states with I ^ 2 of 
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the antiprotonic helium ion can be considered as hydrogen- like, however ns and np-states 
have considerable complex energy shifts. Dependence of AE n i on n is well known from the 
theory, AE ns = AEi s /n 3 , AE np = 32(n 2 — l)/(3n 5 ) • AE 2p . The remaining parameters are 
taken as a combination of the experimental data and theoretical calculations using optical 
model. For 4 He we use Ti s = 11 keV, T 2p = 36 eV, e n i = 0.3r n ;. 
Let us introduce the basis function 

$f /w (r, n B ) = (<Mr) ® Y L (n R )) JM , (10) 

where L is the quantum number of the relative angular momentum of the colliding sub- 
systems, 7T = (— l) l+L is the parity of the system, and J, M are the quantum numbers of 
the total angular momentum (J = 1 + L). Multi-index k involves the set of the quantum 
numbers {n,l,L}. In addition, we introduce also the functions Q{ M7T (r, Qr), which con- 
tain R n i(r) and form with ({TO]) a bi-orthogonal system. Total wave function of the system 
(j»He 2+ ) + He in the state with the total energy E, the total angular momentum J and the 
parity tt in the framework of the close-coupling approximation can be written as 

tff JM ^(R, r) = J2 ^(r, n R )i/£ J *{R)/R, (11) 

k 

where i is the index of the incoming channel. Substituting (|lip into the Schroedinger 
equation with the effective hamiltonian (j7j) and taking into account the bi-orthogonality 
property of the basis states, we get the system of coupled-channel equations 

^(R) + [k] - Lj(Lj + 1)/R 2 ] ^(R) = 2mY,V jk (R)ip ki (R), (12) 

k 

where kj = 2m(E — E n {) and Vjk(R) = (<&j\V(R, r)|$fc). For the channels with non-zero 
T n i we assume Imkj > 0. The conserved quantum numbers (E, J, M, n) are omitted the for 
the sake of brevity. 

The system (1121) looks very similar to the commonly used in the standard close-coupling 
methods. However in our problem the wave numbers kj in some channels contain imaginary 
parts, therefore the conventional boundary conditions at R —>■ oo (incoming + outgoing 
waves) in these channels have no sense and should be modified. (Of course, the boundary 
conditions at the origin ^(0) = remain valid.) For simplicity, we will not consider the 
closed channels with E^ > E. Let a be the subspace of the normal channels (with T n i = 0) 
and (3 the subspace of the annihilating channels with T n i ^ and \mkj > 0. The channels 
from the subspace (3 can not be presented among the incoming channels, therefore the total 
wave function \l/f JM7r (R, r) at i 6 (3 has to be zero everywhere in R e [0, oo), i.e., 

ipji(R) = at i G P, any j. (13) 

For the ordinary incoming channels (i G a) the components ijjji(R) with j G a can contain 
the both incoming (at j = i) and outgoing waves, whereas the components with j G /3 have 
to damp out at the large distance. They also can be presented in the form of 'outgoing' 
waves with the complex kj (Imkj ^ 0). Hence the asymptotic boundary conditions at 
R — > oo for the functions tpji(R) at j 6 a and any j can be written as 

ipji(R) -> k; 1/2 exp[-i(kiR - L l n/2)}5 ij - kj 1/2 exp[i(kjR - L j n/2)}C ji at % G a, (14) 
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For the channels j G a this asymptotic expression has the traditional sense, and the coeffi- 
cients at the outgoing waves give the elements of S-matrix, Sji = Cji, corresponding to the 
transitions between the channels in the subspace a. For the channels j G ft the first term 
in (1T2T) is absent because of j ^ i, whereas the second term is, in fact, damped out: 

ipji(R) ~ exp(—lmkj ■ R)Cji at i G a, j G ft. (15) 

The density current of the outgoing waves in these channels is also damped out exponen- 
tially, therefore the coefficients at j G ft, i G a can not be treated as the S'-matrix 
elements of the transitions (a — > ft). 

The S'-matrix in the subspace a gives elastic and inelastic cross sections 

a(nl -> n'l') = ^ ^ 2^+1 \( n ' l ' L '\ sJ \ nlL ) ~ ^nn'h'^LL'f • (16) 

* L,L',J 

Due to absorptive part of the hamiltonian h the S-matrix is non-unitary The defect of 
unitary gives the cross section of induced annihilation for the definite initial state nl, 



7T 2J+ 1 



k?j-j(2l + l)(2L + l] 



1 - \(n'l'L'\S J \nlL)\< 



n'l'L' 



(17) 



More detailed discussion of the problem of boundary conditions for the close-coupling 
method with annihilating states included in the basis, as well as the scheme of the numerical 
solution of the system (Tl2|) satisfying to the boundary conditions ^(0) = at the origin, 
to the constraint (fl3|) at i G ft and to the boundary conditions (Tbfl) at R — > oo, will be 
given in the separate paper [7]. 



4 Stark transitions and induced annihilation in an- 
tiprotonic helium ions 

With interaction ([3]) and complex energy shifts of ns and np-states we have solved the 
coupled-channel equations using the basis that includes all states with different I at fixed 
n and relative angular momenta up to L = 12 at E ~ 10 K. For the calculation of matrix 
elements Vjk(R) we use the additional approximation R n i(r) = R n i(r) = Rni~ hke ( r ) for all 
states included the basis. 

Cross sections of elastic scattering, of partial and total Stark transitions and of induced 
annihilation for the initial states with / ^ 2 were obtained. A relative importance of the 
different terms in the interaction can be illustrated by the total cross sections of the Stark 
transitions from the circular state with n = 30 at E = 10 K. If we take the short-range 
Morse potential for Vo(-R), the cross section is about 120 a.u. When the polarization term 
V P (R) is included into Vo(-R), the same cross section rises to 620 a.u., whereas the quadruple 
term in ([3]) changes the cross section by the negligible value (< 0.2%). Therefore the main 
part of the calculations was done without the quadruple interaction. 

Figure [T] shows the dependence of total Stark and induced annihilation cross sections on 
the initial quantum number / of the antiprotonic ion at n = 30 and E = 10 K. The induced 
annihilation has an appreciable cross sections for / < 15. At higher initial /, this process is 
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Figure 1: Dependence of total Stark 
(curves 1 and 2) and induced annihi- 
lation (curve 3) cross sections on ini- 
tial I (n = 30, E = 10 K). The curve 1 
is obtained without including the an- 
nihilation widths. 

negligible and doesn't affect the Stark transitions. It should be noted that the total Stark 
cross sections at these energies are very large, of several hundreds atomic units. They are 
comparable with the cross section (j2J) of 'polarization capture' (a p ~ 930 a.u. at E = 10 
K). Fig. [2] shows the energy dependence of the total Stark cross sections for circular orbits 
with n = 30 for two isotopes 4 He and 3 He. The curves look similar, but don't give simple 
isotopic scaling. Dependencies of the Stark cross sections on the principal quantum number 
n for two isotopes at the cm. energy E — 10 K are illustrated by Table [H 

Table 1: Total Stark cross sections (in a.u.) for the circular orbits at cm. energy E = 10 
K depending on the principal quantum number n for two isotopes. 





n 


isotope 


28 


29 


30 


31 


32 


4 He 


573.5 


598.2 


619.6 


636.5 


649.0 


3 He 


539.3 


566.5 


583.4 


597.0 


605.1 



Figure 2: Energy dependence of the Stark 
cross section for the circular orbit with n = 
30 for 4 He (full curve) and 3 He (dashed 
curve). 



Partial Stark cross section a{nli — > nlf) for 4 He, n = 30, at E = 10 K are shown on 
Fig. (Elastic cross sections Al = are not shown here, because they are out of the 
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Figure 3: Partial Stark cross sections a(nli — > nlf) for collisions (p 4 He 2+ ) n ;+ 4 He, at n 
E = 10 K. 



30, 



boundary of the figure.) It is seen from the figure that the cross sections tend to diminish 
with Al = \li — lf\. Nevertheless, in spite of the fact that the coupling between channels 
is carried out by the dipole interaction, the cross sections remain appreciable up to Al of 
several units due to multi-step transitions that are included in the close-coupling equations. 



Table 2: Per-atom Stark transition rates at T = 10 K for circular orbits of two helium 
isotopes (units 10~ 16 MHz • cm 3 ). 



isotope 


n = 28 


n = 30 


n = 32 


4 He 


3.77 


4.15 


4.46 


3 He 


4.46 


4.82 


5.01 



Per-atom rates of Stark transitions from circular orbits, Kst = (vest), averaged over the 
thermal motion at T = 10 K, are given in Table [2] for the two isotopes. The calculated 
values can not be compared directly with the experimental values [3] of per-atom effective 
annihilation rates d^/dp ~ (1 — 3) • 10~ 16 MHz • cm 3 , because the latter rates are a result 
of many cascade processes. Nevertheless, Stark transitions from circular orbits can be 
considered as the most important first step of the cascade, because they change the inner 
angular momentum of the antiprotonic ion and open up possibilities for other processes. It 
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is important that the values of ust are greater than (but of the same order of value as) the 
experimental rates dj/dp. Moreover, we obtain also that K St for 3 He is greater than for 
4 He and grows with n, and these tendencies are correlated with the observed dependencies 
of the value d^/dp. 

5 Conclusion 

We have considered Stark transitions and induced annihilation in the collisions of cold 
(joHe 2+ ) n ; ions with He atoms, using the coupled-channels method with the annihilating 
states (ns and np) included in the basis and with the respectively modified boundary 
conditions. The results of the calculations for the states with n ~ 30 show that the induced 
annihilation is appreciable for the initial states with orbital angular momenta I < 10 — 15. 
For the larger /, the induced annihilation is negligible and the Stark cross sections are not 
affected by the annihilating states. Total Stark cross sections are comparable by order of 
value with the classical polarization capture cross section. Partial Stark cross sections are 
appreciable for the transitions with Al up to several units. 

Per-atom rates of Stark transitions K S t from circular orbits for two helium isotopes, 
averaged over the thermal motion, reveal a some increase of the rates with n and the 
isotope effect, which show the (12 - 18)% excess in the rates for 3 He. These effects are 
correlated with the observed dependencies of the effective annihilation rates d^jdp for the 
initial circular orbits. The calculated and measured values, Kst and d^f/dp, respectively, are 
different but closely coupled in the physical content. We have obtained also that the values 
of K St are greater than, but of the same order of value as the experimental rates d^/dp. 
Therefore we conclude that our results are compatible with the experimental data on the 
collisional properties of cold antiprotonic helium ion in the medium. 
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